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Three Spectacular HII-buried-AGN Galaxies from SDSS 
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ABSTRACT 

We present our analysis of the three HII-buried-AGN: SDSS J091053+333008, 
SDSS J121837+091324, and SDSS J153002-020415, by studying their optical 
spectra extracted from SDSS. The location in the BPT diagnostic diagrams of 
the three galaxies indicates that the narrow emission lines are mainly exited from 
HII regions. However, after the removal of the host galaxy's stellar emission, the 
emission lines display the typical feature of Narrow-line Seyfert 1-like. All of the 
three objects have large Eddington ratio, small black hole mass, and low star 
formation rate. We propose that the three galaxies are at the transit stage from 
the starburst-dominated phase to AGN-dominated phase. 

Subject headings: galaxies: active — quasars: emission lines — quasars 

1. Introduction 

In recent years, there have been remarkable evidences that the evolution of supermassive 
black hole (SMBH) and its host galaxy are tightly linked. The significant correlation between 
the mass of the SMBH and the velocity dispersion of the bulge where the SMBH resides in 
implies that there exists a tight relation between AGN activity and star formation in the 
bulge (e.g. Magorrian et al. 1998; Shields et al. 2003; Greene & Ho 2006; Tremaine et 
al. 2002; Ferrarese et al. 2006). It is now clear that AGN activity and star formation 
frequently happen together (e.g. Gonzalez Delgado 2002). A specific example is the "Q+A" 
galaxies, characterized by composite spectra displaying broad emission lines as well as Balmer 
absorption lines of A- type stars (e.g., Brotherton et al. 1999; Zhou et al. 2005; Wang & 
Wei, 2006; Goto 2006; Wild et al. 2007). 
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Accumulating studies show a likely evolution sequence of co-evolution of AGN and star 
formation. Using the large SDSS spectra database, Heckman et al. (2004) found that most 
accretion occurs onto the black holes with high stellar surface mass densities and young stellar 
populations. Wang et al. (2006) found that the well-documented Eigenvector 1 (El) space 
(Boroson & Green, 1992; hereafter BG92) could be extended to the infrared color a(60,25). 
They argued that AGNs might evolve along with the El space. Wang & Wei (2008) recently 
studied a sample of partially obscured AGNs from SDSS. The broad Ha emission inferred 
Eddington ratio (L/L Edd ) is found to decrease with age of stellar population. This result is 
consistent with the studies of Kewley et al. (2006) and Wild et al. (2007), who estimated 
L/L Edd indirectly through L([OIII])/cr* for obscured AGNs. 

However, we still know very little about the details of the co-evolution of the starbursts 
and SMBH. For instance, do AGN activity trigger starbursts (e.g. Goncalves et al. 1999), 
or starbursts trigger AGN activity (e.g. Weedman 1983), or they happen at the same time 
and co-evolve together? There are some suggestions from the works on the highest redshift 
quasar that the black holes form prior to the assembly of the stellar bulges (Walter et al. 
2004; Riechers et al. 2008). Numerical simulations of galaxies merger including SMBHs 
developed a theoretical light curve of central AGN activity (Di Matteo et al. 2005; Springel 
et al. 2005). In this model, for most of the duration of the starburst, the black hole is 
"buried", being heavily obscured by surrounding gas and dust, especially in UV/optical 
bands (Hopkins et al. 2005a,b). If so, we would expect to identify some AGNs with their 
narrow emission lines dominantly exited from H II region, indicating that the central AGNs 
are surrounded by the gas or dust of the H II region. 

In this paper, we report the identification of three "HII-buried-AGN" from SDSS Data 
Release 4 (DR4, Adelman-McCarthy et al. 2006) MPA/JHU catalog^ They are SDSS J091053 
+333008, SDSS J121837+091324, and SDSS J153002-020415. Their redshifts are 0.116±0.001, 
0.078±0.001, and 0.051±0.001, respectively. All of the three objects show typical AGN fea- 
tures, such as power-law continuum, broad emission line of Ha and H/5, and even strong 
optical Fe II emissions. However, by fitting the narrow emission lines, we can determine 
that all of the three objects are located in the H II regions in the BPT diagnostic diagrams 
(Baldwin et al. 1981). The paper is structured as follows. We present the objects selection 
in §2, and introduce the spectral analysis in §3. The results and discussions are presented in 
Section 4 and 5, respectively. The A cold dark matter (ACDM) cosmology with Q m = 0.3, 

= 0.7, and H = 70 km s -1 Mpc -1 (Spergel et al. 2003) is assumed throughout this 
paper. 



The catalogs could be downloaded from 



http://www/mpa-garching. mpg.de/SDSS/ 
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2. Objects selection 

Because the blueshifted [OIII] emission is thought to be an indicator of outflows in 
AGN (e.g. Boroson 2005; Komossa et al. 2008), our primary motivation is to compare the 
different properties of blueshifted [OIII] profile between star-forming galaxies and Seyferts + 
LINERs, and results will be presented in subsequent papers. We start our objects selection 
from SDSS DR4 MPA/JHU catalogues (Kauffmann et al. 2003; Heckman et al. 2004). 
The catalogs contain a set of physical properties of 567,486 emission-line galaxies, including 
Seyferts, LINERS, and star-forming galaxies. To ensure our sample have a high quality in 
spectra, we need our sample have a median S/N larger than 20 across the whole spectrum. 
Especially, we need the [OIII] and FL3 emission lines have S/N larger than 40. About 3200 
star-forming galaxies have been selected according to the criteria. Then we focus on the 200 
star-forming galaxies which show the largest [OIII] blueshift (~100 km s -1 ), and check their 
spectra by eyes carefully. Only three objects among them, SDSS J091053+333008, SDSS 
J121837+091324, and SDSS J153002-020415, display typical AGN features, such as broad 
Balmer emission lines and Fell complex, superposed on spectra of young stellar population 
stars. The evident AGN features motivate us to carefully re-examine the spectra properties 
for these three objects. 



3. Spectral analysis 

The spectral reduction is performed by standard IRAF procedures. By assuming the 
Galactic extinction curve with R v = 3.1, each spectrum is corrected for Galactic extinction 
correction following the reddening maps of Schlegel et al. (1998). Each of the extinction- 
corrected spectrum is transformed to the rest frame, along with /c-correction, according to 
the corresponding redshift provided by the SDSS pipelines. The spectra at the rest frame 
are displayed in Figure 1. 

For SDSS J091053+333008, because its star light is so weak that the spectrum is dom- 
inated by AGN's featureless continuum and broad emission lines, we fit the continuum by 
the combination of a free powerlaw and the Fe II template provided by BG92. The FWHM 
of the Fell template is fixed to that of its H/3 broad component (e.g., BG92; Hu et al. 2008). 
For the other two objects, the star light components are removed by the principal component 
analysis (PC A) technique (e.g. Li et al. 2005; Hao et al. 2005; Wang & Wei 2008). Briefly, 
a library of pure stellar absorption-line spectra (eigen-spectra) is built by applying the PCA 
method on the single stellar population (SSP) models developed by Bruzual & Chariot (2003, 
hereafter BC03). The range of ages for these models is between lx 10 5 and 2x 10 10 yr. In our 
fitting progress, several components are adopted and listed as follows: the first seven eigen- 
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spectra, an Fell complex template from BG92, and a Galactic extinction curve (Cardelli 
et al. 1989). A \ 2 minimizing algorithm is performed to determine the stellar absorption 
spectrum for each galaxy over the wavelength range from 3700-6800 A, except for the re- 
gions around the strong emission lines: Ha, H/3, [NII]AA6548, 6584, [OIII]AA4959, 5007, and 
[OII]A3727. The left column of Figure 1 shows the results of the AGN continuum/starlight 
removal, the modeled Fell complex, and the isolated AGN emission-line spectra. 

The AGN emission-line spectra are then modeled by the multi-components modeling 
procedure (SPECFIT task in IRAF package) as introduced by Kriss (1994). Each emission 
line is modeled by a set of Gaussian profiles. We model each of the Ha, H/3 and [OIII] emis- 
sion line profiles by two Gaussians: one broad and one narrow components. The width of 
[OIII]A5007 n is set to be the same as that of H/3„. Both intensity ratios [O III]A5007/A4959 
and [N II]A6583/A6548 are fixed to be 3.0. The profile modelings for the Ha and H/3 regions 
are illustrated in Figure 1 (middle and right columns). The results of the line profile mod- 
elings are shown in Table 1. All the uncertainties given in Table 1 are caused by the profile 
modelings. 



4. RESULTS 

4.1. Star-forming origin for narrow emission lines 

Basing upon a set of line ratios, such as [OIII] /H/3, [Nil] /Ha, [SII]/Ha, and [OI] /Ha, 
the BPT diagrams are commonly used as a powerful tool to determine the dominant energy 
source in the emission-line galaxies. Comparing with star-forming galaxies, AGNs (Seyferts 
and LINERs) have larger line ratios because of their stronger and harder ionizing field. 
The theoretical demarcation lines determine the lower limit of AGNs (Seyfert and LINERs) 
were proposed by Kewley et al. (2001) according to their theoretical stellar photoionization 
models. Using the large spectra database provided by SDSS, Kauffmann et al. (2003) 
proposed an empirical demarcation line separating star-forming galaxies in the [OIII] /H/3 
vs. [Nil] /Ha diagram. 

Because the Balmer line profiles are modeled by a combination of a broad and a narrow 
components in the three objects, it is therefore necessary to re-examine their locations on 
the BPT diagrams in the current study. The three objects are plotted in the BPT diagrams 
in Figure 2. The four panels represent the four different pairs of line ratios. It is clear that 
all the three objects are below the theoretical demarcation lines (see Panel A, B and C) 
proposed in Kewley et al. (2001). In particular, SDSS J091053+333008 is obviously located 
in the HII region in the [O III] /H/3 vs. [N II] /Ha diagram, i.e, far below the empirical 
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classification line suggested by Kauffmann et al. (2003). In the same diagram, the other two 
objects are located marginally above the empirical demarcation line. Since the objects are 
extracted from star- forming galaxies from the MPA/JHU catalog, the reason why they are 
not located in the HII region is that we have removed the broad component of Ha and H/5. 
As a result, the ratio of [OIII]/H/3 and [Nil] /Ha might be a little different from Kauffmann 
et al. (2003). 

As an additional test, we plot the three objects on the [O III]/[0 II] vs. [O I]/Ha 
diagram in Panel D in Figure 2. The line ratio of [O III]/[0 II] is sensitive to the ionization 
parameter of the gas (Kewley et al. 2006). To obtain the line ratio, we need to correct the 
intrinsic dust extinction on the flux of [O III] and [O II]. For each object, the extinction 
is estimated from the observed Balmer decrement for narrow components of Ha and H/5, 
assuming the Case B recombination. The predicted Ha/H/3 ratio is 2.86 for H II region in the 
condition of gas density 100cm -3 and temperature 10 4 K (Osterbrock 1989). The derived A v 
values [=R„xE(B-V)] are listed in Table 1 for the three objects. As shown in the diagram, 
all of our three objects are located in the region occupied by star-forming galaxies. 

In summary, using the widely used diagnostic tools, we draw a conclusion that the 
narrow emission lines in the three objects are mainly powered by the ionizing field from hot 
stars rather than AGNs. 

After demonstrating that the narrow emission lines are mainly powered by star formation 
activity, we attempt to estimate current star formation rate for the three objects in terms 
of the [O II]A3727 emission. [O II]A3727 is an empirical indicator to estimate the current 
star formation rate (SFR) for emission-line galaxies (e.g. Gallagher et al. 1989; Kennicutt 
1998; Kewley et al. 2004). We calculate the SFR according to the equations in Kewley et 
al. (2004): 

SFRrtOTm - 7.9 x 10- 42 L([OII])(ergss- 1 ) 

SFR([0II]) " (-1.75)[log(0/H) + 12] + 16.73 M@yi (1) 
The metallicity is estimated according to the empirical calibration using the R23 indicator 
(Kewley et al. 2004), where R 23 = log([OII]A3727 + [OIII]A5007)/H/3. Again the [O II] 
luminosity is corrected for the intrinsic extinction derived from the Balmer decrement. The 
calculated current SFRs of the three objects are listed in Table 1. It turns out that all of 
them have small SFR~ 1M yr -1 , which indicates a suppressed star formation activity likely 
due to the AGN feed back (see below for discussion of AGN properties). 
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4.2. 



Broad emission lines: an indicative of AGNs 



The AGN properties are examined in this section. The detection of broad Balmer 
emission lines is an indicative of the existence of AGN in the center of galaxy (see Section 4 
for the discussion of other origins). Our spectral measurements indicate that all of the three 
objects show AGN line emission with relatively narrow H/3 profiles and large Fe II ratios, 
which are the typical features of Narrow-line Seyfert 1 Galaxies (NLSls, e.g., Osterbrock & 
Pogge 1985). The measured FWHMs of H/3 broad emission lines are around 2000km s _1 , 
and the Fe II ratios RFe ~ 2 (except SDSS 153002-020415), where the RFe is defined as the 
flux ratio of optical Fe II/H/Sb- The flux of the Fe II blends is measured between the rest 
frame wavelength 4434 and 4684A. Especially, the object SDSS J121837+091324 has already 
been included in the ~2000 NLS1 galaxies sample from SDSS by Zhou et al. (2006). 

The co-evolution of AGN and its host galaxy implies that young AGNs are associated 
with young stellar population and intensive star formation activity. Mathur (2000) proposed 
that NLSls with small black hole mass (Mbh) and high Eddington ratio (L/L-^^) are at 
early stage of evolution of AGN. The evident broad emission lines allow us to study the black 
hole accretion directly in the three objects. The black hole mass (Mbh) and Eddington ratio 
(L/L E dd) are two basic physical parameters of AGNs (e.g. BG92; Boroson 2002). We refer 
the readers to McGill et al. (2008) for a summary of the existing formula used to estimate 
M B h basing upon "single-epoch" spectrum. These methods are based on the great progress 
in the calibration of the -Rblr — L relationship (e.g., Kaspi et al. 2007 and references therein). 
In the current study, we estimate the Mbh for the three objects according to their broad Ha 
and H/9 emission lines (Greene & Ho 2005): 



Lboi/LEdd is estimated by simply assuming the bolometric luminosity of quasar is pro- 
portional to the luminosity at rest-frame wavelength 5100A: L bol = 9AL A (5100A) (e.g., Kaspi 
et al. 2000), where 



For each object, both the luminosities of Ha and H/3 are corrected for local extinction 
according to their observed Balmer decrements. The estimated M B h from Ha and H/5 are 
listed in Table 1, and they are highly consistent with each other within a factor of two 






(3) 
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for all the three galaxies, which is less than the typical absolute uncertainties of the mass 
scaling relationships (Vestergaard & Peterson 2006). It is clear that all the three objects 
show low M B h ~ 10 6 M Q and high L/L Edd ~ 0.2, which implies that the three objects have 
an NLSl-like nucleus in their galaxy centers. 

5. Discussion 

The spectral analysis given above clearly reveals that the emission lines of the three 
objects consist of two components with distinct origins. The narrow emission lines are 
mainly powered by star formation activity, and the broad emission lines by central AGN. 
In addition, the AGNs show small M B h, large L/L Edd , and strong Fe II emission, which is 
typical of NLS1 galaxies (e.g., Zhou et al. 2006). 

Except for the AGN's contribution, other mechanisms may be responsible for the broad 
line emission in the spectra of these objects, including stellar winds and supernova remnants 
(e.g. Izotov et al. 2007). Both of these mechanisms can produce broad emission lines 
with FWHM > 1000 km s -1 . However, the luminosity of Ha broad components for our 
three objects are ~ 10 42 erg s" 1 , far exceeding the typical luminosity of broad Ha of stellar 
winds from Wolf-Rayet or luminous blue variable stars (10 36 -10 40 erg s _1 ). Moreover, if the 
broad component is caused by supernova events, the magnitude obtained from the spectrum 
is expected to greatly differ from the photometric magnitude, because the brightness of 
supernova must decay at the time interval between the two observations of about a year. 
We obtain the SDSS r' band magnitudes from the observed spectra by convolving each 
spectrum with the r' band transmission curve, and compare them with the corresponding 
fiber magnitude obtained in the SDSS The magnitude difference is found to be less than 1 
mag within the interval of about a year for these three objects. As a result, we can eliminate 
the possibility of Type Hp supernova events. 

We propose that the three objects are at the transition stage from the starburst- 
dominated phase to AGN-dominated phase, which is supported by the theoretical evolu- 
tionary scenario given by the numerical simulations of mergers of two gas rich galaxies. 
Hopkins et al. (2005a, 2005b, and references therein) proposed a scenario that the central 
black hole of AGN is buried by gas or dust in the young AGN phase. In their analysis, when 
starbursts are triggered by the nuclear inflow and feed black hole growth, the black hole is 
heavily obscured by surrounding gas or dust. As the black hole grows, the black hole feed- 
back will blow the gas or dust away in a powerful wind, and the luminous AGN will be seen 
at that time. A case for the model is NGC 6240, of which the X-ray data reveal that there 
are two AGNs, while the optical data cannot resolved them, possibly for the reason that the 
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AGNs are in highly dust-enshrouded environments (e.g. Komossa et al. 2003; Gressen et 
al. 2004; Max et al. 2005). We believe that the next generation hard X-ray missions with 
enhanced capability (in sensitivity and imaging), such as NuSTAR, Simbol-X and NeXT, 
would be helpful in identifying and studying the starburst obscured AGNs. 

Finally, we need to understand why the emission from AGN NLR is so weak that 
the narrow emission is overwhelmed by the contribution from H II region. Netzer et al. 
(2004, and see also in Netzer et al. 2006) proposed that the NLR would be dynamically 
unbounded if the AGN luminosity is high enough. The distance of the NRL, if it exists, 
from the central black hole is estimated to be not far less than lkpc for the three objects 
according to the relationship -Rnlr = 1-15L^ 9 42 kpc (Netzer et al. 2004), where Lh/3,42 = 
Z/H/3/IO 42 erg s _1 . At this distance, the escape velocity for a spherical galaxy roughly equals 
to 290M*( 2 1 i?~ 1 k / p 2 c km s _1 , where M*. n is the stellar mass in units of 10 11 M o , and -Riokpc 
the radial distance in units of lOkpc. The stellar mass could be inferred from the the 
estimated mass of the SMBH according to the well-established black hole and bulge mass 
relation (Mbk/M*) ~ 0.002 (McLure & Dunlop 2004), which results in an escape velocity 
only ~ 100 km s -1 . This scenario is supported by recent observation study of extended 
emission-line regions (EELRs) around low-z QSOs. Husemann et al. (2008) found that the 
presence of EELRs are preferentially detected in QSOs with large black hole masses. 

The lack of typical NLR in the three objects leads us to suspect that the NLR might 
evolve with AGN and its host galaxies. Netzer et al. (2004) proposed a scenario of "star- 
forming NLRs" in which the compact NLR gas in kpc scale is caused by the star formation 
activity. An alternative possible scenario proposed here is that the the formation of NLR is 
the consequence of the feedback of central AGNs. As discussed above, a powerful feedback 
wind is required for AGNs to dispel the surrounding gas and dust material to appear as a 
luminous AGN at about 1 Gyr since the merger (Hopkins et al. 2005a). This implies the 
typical NLR would not emerge until a majority of surrounding material has been dispelled. 
Another implication for the wind scenario is that the NLR kinematics is dominated by radial 
outflow. In fact, spectroscopic studies with high spatial resolution of NLR of a few nearby 
Seyfert 2 galaxies indicate that radial velocity fields can be modeled by biconical outflows 
where the NLR gas accelerates outward from the nucleus to a turnover point at a distance 
of ~100pc from the nucleus and subsequently appears to decelerate to the systemic velocity 
at a distance of ~300pc (e.g., Hutchings et al. 1998; Kaiser et al. 2000; Nelson et al. 2000; 
Crenshaw et al. 2000; Crenshaw & Kraemer 2000; Ruiz et al. 2001). 
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6. Conclusions 

By studying the optical spectra of the three objects: SDSS J091053+333008, SDSS 
J121837+091324, and SDSS J153002-020415, we find that all of them are located in the HII 
region in the BPT diagnostic diagram. However, all of their spectra display characteristic 
features of AGNs. All of them have large Eddington ratio, small black hole mass, and small 
SFR. We propose that the dust-buried-AGN (Hopkins et al. 2005a, b) can explain the 
feature. 

We would thank the anonymous referee for constructive comments that improved the 
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template. This work was supported by the National Science Foundation of China (grant 
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Fig. 1. — left: Detailed modeling and subtractions of the starlight components in the spectra 
of the three objects. For SDSS J091053+333008, we plot the observed spectrum, emission- 
line spectrum, a power-law continuum, and Fell blends, from top to bottom. For the other 
two objects, we plot the emission-line spectrum, observed spectrum, modeled starlight spec- 
trum, and Fell blends. The spectra are vertically shifted by arbitrary amounts for visibility. 
Middle: Modelings of the two Gaussion profiles for the H/5 region. Right: Modelings of a set 
of Gaussion profiles for the Ha region. 
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Fig. 2. — The four BPT diagnostic diagrams for our three galaxies. The theoretical lines 
separating AGNs from star-forming galaxies proposed by Kewley et al. (2001) are shown by 
the solid lines in all four panels, while the empirical lines proposed by Kauffmann et al. (2003) 
by the dashed line in Panel A. The dot-dashed lines in Panel B, C, D separating LINERS from 
Seyferts are proposed by Kewley et al. (2006). The circles are for SDSS J091053+333008, 
the squares are for SDSS J121837+091324, and the triangles are for SDSS J153002-020415. 
All of our three galaxies are located in the star-forming region, which is proposed by Kewley 
et al. (2006). 
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Table 1: Properties of the emission lines of the three objects, Black hole masses, Eddington 
Ratios and SFR. 



Lines 


SDSS J091053+333008 


SDSS J121837+09 


1324 SDSS J153002-020415 


H/3 N 


25.8±1.3 


11.1±0.7 


27.5±2.0 




37.3±1.7 


14.8±0.9 


39.5±2.8 


[OIII]A5007 


20.3±0.7 


16.1±0.4 


41.1±0.5 




99.5±3.1 


45.5±1.5 


93.0±4.1 




148.1±4.0 


79.0±2.2 


185.8±6.4 


[NII]A6583 


25.3±1.2 


15.9±0.7 


30.2±2.3 


[SII]A6716 


10.0±0.6 


7.8±0.4 


17.4±1.4 


[SII]A6731 


9.5±0.6 


5.2±0.4 


11.8±1.2 


[OI]A6300 


2.2±0.5 


2.4±0.5 


4.1±0.9 


[OII]A3727 


17.4±1.0 


14.8±0.7 


65.5±3.5 


Rfe 


1.7±0.1 


2.9±0.2 


O.liO.O 


FWHM(H/3 B ) 


2.3±0.2 


1.7±0.2 


1.7±0.2 


A v 


0.70±0.04 


0.88±0.04 


0.35±0.02 


M BH (Ha)/10 6 M o 


7.9±2.0 


3.2±0.7 


2.8±0.6 


M BH (H/3)/lO 6 M 


9.9±2.4 


2.3±0.7 


1.7±0.6 


Lbol/LEdd 


0.21±0.04 


0.16±0.03 


0.13±0.03 


SFR 


2.4±0.1 


0.8±0.0 


0.6±0.0 



Notes: The flux of each component is in the unit of 10 16 ergs s 1 cm 2 . The FWHM(H/9) 
is in the unit of 10 3 km s -1 . And the SFR is in the unit of M yr _1 



